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Abstract 

Autocatalysis is a fundamental concept, used in a 
wide range of domains. From the most general def- 
inition of autocatalysis, that is a process in which a 
chemical compound is able to catalyze its own for- 
mation, several different systems can be described. 
We detail the different categories of autocatalyses, 
and compare them on the basis of their mechanistic, 
kinetic, and djmamic properties. It is shown how 
autocatalytic patterns can be generated by different 
systems of chemical reactions. The notion of auto- 
catalysis covering a large variety of mechanistic re- 
alisations with very similar behaviors, it is proposed 
that the key signature of autocatalysis is its kinetic 
pattern expressed in a mathematical form. 

Keywords chemical network, autoinduction, tem- 
plate, competition, mechanism 
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1 Introduction 

The notion of "autocatalysis" was introduced by Ost- 
wald in 1890^ for describing reactions showing a rate 
acceleration as a function of time. ^ It is for example the 
case of esters hydrolysis, that is at the same time acid 
catalyzed and producing an organic acid.^ Defined as a 
chemical reaction that is catalyzed by its own products, 
it has quickly been described on the basis of a character- 
istic differential equation.^ Typically used to describe 
complex behaviors of chemical systems, like oscillatory 
patterns,^ it has immediately appeared to be essential for 
the description of biological systems: growth of individual 
living beings,^ population evolution^ or gene evolution.!^ 

Extending this concept from a chemical description to a 
more open context was initially carefully described as an 
analogy, sometime qualified by the more general notion of 
"autocatakinesis".' ^° However, this eventually leads to an 
overgeneralization of the term of autocatalysis, tending 
to be assimilated to the notion of "positive feedback", for 
example in economy. 

The notion of autocatalysis is now actively being used 
for describing self-organizing systems, namely in the field 
of emergence of life and artificial life. Autocatalytic pro- 
cesses are the core of the mechanisms leading to the sym- 
metr y brea king of chemical compounds towards homochi- 
rality^ Fi nland could be identified in several experimental 
systems IHEE However, how such autocat alytic processes 
shall manifest is still under heavy debateli^iiZI 

The purpose of this article is thus to clarify the meaning 
of chemical autocatalysis and this effort will be under- 
taken by covering these following points: 

• What is autocatalysis for a chemical system? On the 
basis of the general description of a process allow- 
ing a chemical compound to enhance the rate of its 
own formation, autocatalysis is defined by a kinetic 
signature, expressed in a mathematical form. 

• How can an autocatalytic process be realized? As 
many mechanisms can reduce to the same macro- 
scopic kinetic laws exhibiting autocatalysis, the focus 
is put on several mechanistic realisations of auto- 
catalytic processes, based on simple models further 
illustrated by concrete chemical examples. 

• How can autocatalysis be observed and character- 
ized? The focus is put on the djTiamic properties, 
showing that this observable is the direct conse- 
quence of the kinetic pattern, rather than the un- 
derljdng mechanism. 
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2 Autocatalysis: a Practical Definition 



What is the role of autocatalysis? Embedded in non- 
equilibrium reaction network, the competition be- 
tween autocatalytic processes allows the onset of 
chemical selection, that is the existence of bifurca- 
tion phenomena allowing the extinction of some com- 
pounds in favor of others. 



A, and is thus more important at the beginning of the re- 
action (thus an initial exponential increase of the product 
B) that at the end (thus a damping of the autocatalysis) 
resulting in a global sigmoidal evolution. In systems were 
the influence of A on fc is weaker, as detailed further, an 
undamped autocatalysis will be observed characterized by 
an exponential variation until the very end (see fig.[2][c)). 



2 Autocatalysis: a Practical 
Definition 

2.1 A Kinetic Signature 

From its origin, the notion of autocatalysis has focused 
on the kinetic pattern of the chemical evolution.^ The 
general definition of autocatalysis as a chemical process 
in which one of the products catalyzes its own formation 
can be mathematically generalized as: 

^ = fc(X)-x;'+/(X), fc>0; n>0; |fc|»|/| (1) 
at 

X is the vector of all the concentrations Xj . An autocatal- 
ysis for the compound exists when the conditions of 
eq. ( ij are fulfilled. The term fc(X) • x" describes the auto- 
catalytic process itself, while /(X) describes the sum of all 
other contributions coming from the rest of the chemical 
system. 

We have an effective practical definition of the con- 
cept of autocatalysis, based on a precise mathematical 
formulation. The causes of this kinetic signature can be 
investigated, searching what mechanism is responsible for 
the autocatalytic term. This leads to the discovery of a se- 
ries of different kinds of autocatalysis processes, and their 
respective effect, describing what observable behavior is 
generated by the autocataljrtic term (see fig.[l]). 

2.2 Potential vs Effective Autocatalysis 

This kinetic definition is purely structural. As a mat- 
ter of fact, a system may contain potentiaZ autocatalysis 
i.e. an autocatalytic core exists in the reaction network. 
However, in the absence of some specific conditions nec- 
essary for this autocatalysis to be effective, the potential 
autocatalysis may be hidden by other kinetic effects, and 
thus not manifests its behavior in practice. 

Possibly, in eq. ([l]), the term /(X) may simply over- 
whelm the autocatalytic process. This is typically the 
case when an autocatalysis is present together with the 
non-catalyzed version of the same reaction, that may not 
be negligible in all conditions. A simple example is a 
system simultaneously containing a direct autocatalysis 
A-\-B — > 2B, concurrent with the non autocatalytic reac- 
tion A — > B. The autocatalytic process follows a bimolec- 
ular kinetics, and will be more efficient in a concentrated 
than in a diluted solution. The dynamic profile of the re- 
action is thus sigmoidal for high initial concentration of A, 
but no more for low initial concentration (see fig.|2];^a-b)). 

It can also be seen that the term fc(X) may vary during 
the reaction process. In a simple autocatalytic process as 
described above, k is proportional to the concentration in 



3 Mechanistic Distinctions 

How can this kinetic pattern be realized? Let us now de- 
tail several types of mechanisms. They can all be reduced, 
in some conditions, to the autocatalysis kinetic pattern of 
eq. ([T]). All of them will be equally defined in the paper as 
autocatalytic, while this status may have been disputed in 
the past on account of the distinct chemical realisations. 
In the following, we emphasize the major mechanistic 
pattern to eventually be reduced to an equivalent kinetic 
autocatalysis, and discuss where their difference comes 
from. 

3.1 Template Autocatalysis 

The simplest autocatalysis is obtained by the X — > 2X 
pattern. It can be represented by: 



A + B ^B+B 



(2) 



The corresponding network is given in fig. [3][a). It can 
further be decomposed through the introduction of an 
intermediate compound C : 



A + B- 



■B + B 



(3) 
(4) 



The corresponding network is given in fig. |3](b) . 

The first mechanism entails the following kinetic evolu- 
tion: 



db _ da 
d7~~ dt 



(5) 
(6) 



This can be expressed as a chemical flux = by 



relying on the Mikulecky formalismP^n^S 
a 



b 



Ti — fci • Ka^b — fc-i ■ 



(7) 
(8) 

(9) 

(10) 



fcj and k_i are the kinetic constant rates of the reaction 
1 in the direct and reverse direction. K^^ and Kg are the 
thermodynamic constant of formation of compounds A 
andB. 

Formally there is a linear flux of transformation of A 
into B, coupled to a circular flux of same intensity from B 
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3 Mechanistic Distinctions 



Mechanistic: 

Source of AC 



nX mX 
Strict AC mechanism 




Autoinductive mechanism Other mechanisms? 



Kinetic: 

Definition of AC 



dt 



i 



k-x2 + f{x) 



i 



Dynamic: 

Observation of AC 



Linear scale: Logaritlimic scale: Inverse scale: 

n>0 convex (rate acceleration) n>l convex (over-exponential) n>2 concave (over-hyperbolic) 

n=0 concave (non-autocatalytic) n<l concave (sub-exponential) n<2 convex (sub-hyperbolic) 





Figure 1 - Classification of the concepts of autocatalysis (AC) depending on their descriptions (mechanistic, kinetic, and dynamic). 
The graphs represents the time evolution of a non-autocatalytic reaction (n = 0, red), and of autocatalytic reactions of order n = 1/2 
(green), 1 (blue), 3/2 (dotted red), 2 (dotted green), and 3 (dotted blue). 




(a) No autocatalysis (b) Autocatalysis 



(c) Undamped 
autocatalysis 



Figure 2 - (a-b): First order autocatalytic process (Tj = 10^ 
M.s~^) in presence of a non-autocatalytic reaction (r2 = 
10"^ M.s~^) of spontaneous transformation of A into B 
(Kj, = \ M, Ks = \0'^ M). (a) Diluted (a„ = 10"^ M). 
(b) Concentrated (a^ = 1 M). (c) Undamped autocataly- 
sis (Indirect autocatalysis, described in fig.^p^), = 0.1 

back to B (see fig.|3][a-b)). In presence of an intermediate 
compound, the equations becomes: 



^2 = ^2(yc-VB^ 



(11) 

(12) 



Under the hypothesis that C is an unstable intermediate, 
(i.e. <C KgjKy^), the variation of C can be neglected 
compared to the variations oiA and B (quasi steady-state 
approximation, hereafter QSSA), so that: 



systems, we are in presence of the perfect kinetic signature 
of an autocatal5rtic system, following a sigmoidal evolution 
(see fig.|4](a)). This equivalence is guaranteed as long as 
the compound C remains unstable. When it is not the 
case, the dimeric intermediate C hardly liberates the final 
compound B, which eventually lea ds to an autocatalytic 
process of order 1/2 rather than IjSMI 

Template autocatalysis requires a direct association be- 
tween the reactants and the products. This is typically the 
case of DNA replication, one double strand molecule giv- 
ing birth to two identical double strand molecules, thanks 
to the very selective association of complementary nu- 
cleotides along each strand. More simple examples can 
be found in some biological mechanisms that requires 
autocatalytic processes, for example for the generation of 
chemical oscillation inducing circadian rhytmicity in cells. 
The system described by Mehra et af^ is based on a non 
equilibrium system of association/dissociation of proteins 
forming a large chemical cycle [C — > AC — > AC* — > 
ABC* — > BC* — > C* — > C], maintained by a flux of 
ATP consumption, one cycle consuming and freeing A and 
B. 122! The oscillations are generated by coupling this chem- 
ical flux to an autocatalytic process of phosphorylation 
obeying to the reaction scheme:?^ A-l- C -I- AC* — > 2AC*. 



^1 =^ V>2 

= V> 



(13) 
(14) 

(15) 



The system is strictly equivalent to the direct autocatalysis, 
with an apparent rate r^rj/Cri -I- Fj). With these two 



3.2 Network Autocatalysis 

The direct mechanism of template autocatalysis is con- 
ceptually the simplest framework. It may actually not 
be the most representative class of autocatalysis, and a 
similar kinetic signature can result from more complex 
reaction networks. 
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3 Mechanistic Distinctions 




(e) Iwamura et a/ system 



(f) Collective 



Figure 3 - Reaction network of different autocatalytic pro- 
cesses of spontaneous transformation of A into B (a-d), of 
A + X into AX (e), and of A; into B; (f). The indicated fluxes 
correspond to what is observed within the QSSA. 



3.2.1 Indirect Autocatalysis: 

The autocatalytic effect can be indirect when reactant 
and products never directly interact: 



^B + E 



E^B 
B^D 



(16) 
(17) 
(18) 
(19) 



There is no direct A/B coupling, nor direct 2B formation, 
but the presence of a dimeric compound C. The network 
decomposition of this system (see fig.[3][c)) implies once 
again a linear flux of transformation of ^4 into B, linked to 
a large cycle of reaction transforming B back to B. This 
system is still reducible to anX — > 2X pattern. 

The QSSA for compounds C, D, E allows to express the 
reaction flux as: 



V ■ 



± + ± + Ya + Yb 



(Vj^Vb - V^) (20) 



The details of the calculations are given in appendix. 

When the terms V^/F^ and Vg/F^ are small compared 
to either F^^ or ^ (i.e. when at least one of the two 
reactions of eq. ([T6]l-eq. dTTll) is kinetically limiting), the 
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Figure 4 - Time evolution of compound concentrations for 
different autocatalytic processes of spontaneous transforma- 
tion of A into B (Kyi=lM and = 100 M) in a logarithmic 
scale for concentrations (a-c), or logarithmic scales for both 
time and concentrations (d). K and concentrations are in M, 
times in s, and F in M.s"\ (a): fig.^b), Tj = 1, r2 = lO"*, 
Kc = 0.01; (b): fig.^c), Tj = r2 = = r4 = 10 (except 
the values indicated on the graph), Kq = Kj^ = = 0.01; 



(c):fig.^d), Y^ = F^ = 100,Kc 
flg.m),F^ = \00,T^ = l. 



1, Ke* = 10; (d): 



system behaves like a simple autocatalytic system, with 
if (xa-h before the reaction completion, with a progres- 
sive damping of the exponential growth as long as A is 
consumed. When the term VJF^ is predominant (i.e 
when the reaction of eq. ( [19] ) is kinetically limiting), the 
flux is oc b: the profile remains exponential up to the 
reaction completion, with no damping due to A consump- 
tion. When the term V^/F-^ is predominant (i.e when 



the reaction of eq. (181 is kinetically limiting), the flux is 
if <xa: the autocatalytic effect is lost (see fig.|4fb)). 

Network autocatalysis is probably the most common 
kind of mechanisms. A typical biochemical example is the 
presence of autocatalysis in glycolysis. In this system, 
there is a net balance following the X — > 2X pattern. ATP 
must be consumed to initiate the degradation of glucose, 
but much more molecules of ATP are produced during 
the whole process. While these systems are effectively 
autocatalytic, there is obviously no possible "templating" 
effect of one molecule of ATP to generate another one. 



3.2.2 Collective Autocatalysis: 

More general systems, reminiscent of the Eigen's hy- 
percycles,^® are responsible of even more indirect auto- 
catalysis. No compound influences its own formation rate, 
but rather influences the formation of other compounds, 
which in turn influence other reactions, in such a way that 
the whole set of compounds collectively catalyzes its own 
formation. 

A simple framework can be built from the association 
of several systems of transformation A^ — > B;, each B; 
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3 Mechanistic Distinctions 



catalyzing the next reaction (see fig.[3][f)): 

Ai+B,_^^Bi+Bi_^ (21) 
(with i = {1, 2, 3, 4} and B4 = Bq) 

There are four independent systems, only connected by 
catal5rtic activities. 

If the system is totally symmetric, then all b,- are equal, 
and all a, are equal, so that the rates become: 



v, = r,VB,_^(VA,-VB,) 
if = rVs(.v^-Vs) 



(22) 
(23) 



This leads to a collective autocatalysis with all compounds 
present. They mutually favor their formation, which re- 
sults in an exponential growth of each compound (see 
fig.|4][d) dotted curve). 

With symmetrical initial conditions (i.e. identical for the 
four systems), the system strictly behaves autocatalytically. 
If the symmetry is broken, e.g. by seeding only one of the 
B;, the system acts with delays. The evolution laws are 
sub-exponential, of increasing order; at the very beginning 
of the reaction, considering that ^4; do not significantly 
change and that B; are in low concentrations, we obtain 
oc t'~^. Seeding with Bj, the compound B2 evolves in 
t^. Its impact on compound B3 induces an evolution in 
t^. In its turn, the impact of compound B3 on compound 
B4 induces an evolution in t'*. The compound 1 at first 
remains constant, and it is only following a given delay 
that it gets catalyzed by B4 (see fig.Qd)). 

This system is actually not characterized by a direct 
cyclic flux, but by a cycle of fluxes influencing each other 
and resulting in a cooperative collective effect: 

(Ai+A2+A3+A4) + [Bi+B2+B:i+B4) (24) 
2(Bi-|-B2-l-B3-|-B4) (25) 

The simultaneous presence of all different compounds is 
needed to observe a first order autocatalytic effect. Given 
asymmetric initial conditions, a transitory evolution of 
lower order is first observed, until the formation of the 
full set of compounds. 

A typical example of collective autocatalysis is observed 
for the replication of viroids.^' Each opposite strand of 
cyclic RNAs can catalyze the formation of the other one, 
leading to the global growth of the viroid RNA in the 
infected cell. 



3.2.3 Template vs Network Autocatalysis: 

All the preceding systems can be reduced to aX — > 2X 
pattern. This is characterized by a linear flux of chemical 
transformations, coupled to an internal loop flux: for each 
molecule (or set of molecules) A transformed into B, one 
B is transformed and goes back to B, following a more or 
less complex pathways. They can be considered as mech- 
anistically equivalent: a seemingly direct autocatalysis 
may really be an indirect autocatalysis once its precise 
mechanism is known, decomposing the global reaction 
into several elementary reactions. 

Practically, autocatalysis will be considered to be direct 
(or template) when a dimeric complex of the product is 



formed (i.e. allowing the "imprint" of the product onto the 
reactant). If such template complex is never formed, we 
preferentially speak of network autocatalysis, in which the 
X — > 2X pattern only results from the reaction balance. 

3.3 Autoinductive Autocatalysis 

Some reactions are not characterized by a X — > 2X 
pattern, but still exhibit a mechanism for the enhancement 
of the reaction rate by the products. This is typically the 
case for systems where the products increase the reactivity 
of the reaction catalyst rather than directly influencing 
their reaction production itself These systems still possess 
the kinetic signature of eq. ([T]), but are sometime referred 
as "autoinductive" instead of "autocatal5rtic".E3 

3.3.1 Simple network: 

Let us take a simple reaction network of a transforma- 
tion A — > B catalyzed by a compound that can exist under 
two forms E/E*, E* being the more stable one. These two 
forms of the catalyst interact differently with the product 
B (see fig.igd)): 



A + E^C 



C^B + E 



C^B + E* 



(26) 
(27) 
(28) 



There is no dimeric compound in the system, even indi- 
rectly formed. 

Provided the catalyst, present in C, E, and E*, is in low 
total concentration, the QSSA implies the presence of two 
fluxes: the transformation of A into B catalyzed by E of 
intensity (f, and the transformation of E* into E catalyzed 
by B of intensity e, with » e. Assuming that E* is very 
stable compared to E and C, this decomposition gives (see 
appendix for details) : 



(29) 



The autoinduction is kinetically equivalent to the indi- 
rect autocatalysis mechanism: 

• When r, » ri§i, the flux w is r-,V^,(V. - Vg): the 
system is non-autocatalytic. 

. When r2 ^ Ti |, the flux >f is Ta^CV^Vg - V^): the 
system is simply autocatalytic. 

. When r2 < , the flux ^ is rjVO (^b - f ) : 
the system presents an undamped autocatalysis. 

Following the kinetic analysis, the behavior is similar to 
the time evolution of autocatalytic systems (See fig.|4](c)). 
The behavioral equivalence of these two systems (kinet- 
ically equivalent but mechanistically very different) will 
be investigated in more details in the next section. 
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4 Embedded Autocatalyses 



3.3.2 Iwamura's model 

An example of autoinductive autocatalysis is the proline- 
catalyzed a-aminoxylation of aldehydes. The core princi- 
ple is a reaction j4-I-X — *AX, catalyzed by P, the product 
AX catalyzing the first catalytic step P +A — > PA (see 
fig.[3][e)). This chemical system can be decomposed into 
two different fluxes A + X — *AX, one coupled to a cat- 
alytic cycle [P — > PA — > PAX — * P\AX — > P], and 
one coupled to a catalytic cycle [PA — > PAX — > P|AX 
— > PA] . The first one contains the slow reaction of A 
on P, and corresponds to a slow flux e. The second one 
only contains fast reactions, and corresponds to a fast flux 
if. In an ideal case (see appendix for details), the flux of 
production of AX is equal to: 



VaVax 



AX 
Vy 



(30) 



The kinetic signature of an undamped autocatalysis is 
once again obtained. 

3.3.3 Network vs Autoinductive Autocatalysis: 

Autoinductive autocatalysis is mechanistically differ- 
ent from network or template autocatalysis. The balance 
equation is rather of the form A+ aB — > (1-1- a)B, with 
a <C 1. The linear transformation A — > B is only weakly 
coupled to the cycle of B back to itself, this latter one being 
subject to a much lower flux than the linear flux. However, 
autoinduction is kinetically and dynamically equivalent 
to network autocatalysis, leading to the same kind of dif- 
ferential equation, and thus of behavior It must be noted 
that the undamped exponential profile — due to a flux only 
proportional to the products and not to the reactant — is 
not characteristic of autoinductive processes^^ but can 
also be explained by network autocatalytic mechanisms, 
when the consumption of the reactant is not limiting the 
kinetic of the network. 



4 Embedded Autocatalyses 

Autocatalysis is not so important per se but as a way 
of giving birth to rich non-linear behaviors like bifurca- 
tion, multistability or chemical oscillations. It is crucial 
to study the interaction of autocatal5rtic mechanisms and 
their ability to generate such behaviors when embedded 
in a larger chemical network. 



4.1 Dynamical Distinctions 

Different behaviors depending on the order n of the 
autocatalysis can be observed in biochemical competitive 
syste ms. T hey are classically studied in population evo- 
lutiori'^°'^^' and described as "survival of the all" in the 
case of < n < 1 (characterized by the coexistence of 
all compounds), as "survival of the fittest" in the case of 
n = 1 (when the only stable solution retains the fittest 
compound or the most "reproductible") and as "survival 
of the first" in the case of n > 1 (when the final solution 



just retains the product initially present in the highest 
concentration) . 

The case < n < 1 is the least interesting one, as it 
hardly leads to a clear selectionnist process. However, 
real mechanism that seems to possess a first order auto- 
catalysis may actually present a lower autocatalytic order. 
This is typically the case for direct template autocatalysis, 
in which the order falls to 1/2 on account of the high 
stability of the dimeric intermediate — ^which is actually a 
necessary condition for the selectivity of template replica- 
jjQjjj2ij|22|32 Yj^jg turns out to be a fundamental problem 

for understanding the emergence of the first replicative 
molecules. I^^HIII 

More complex mechanisms may lead to higher orders, 
t)^ically by the formation of dimeric autocatalysts.^^ This 
is the case of the Soai reaction whose high sensitivity 
to initial conditions may potentially be explained by the 
formation of trimeric^^ or even hexameric complexes. 122! 

4.2 Comparative Efficiency of Direct 
and Autoinductive Autocatalyses 

The relative efficiency of two different autocatalytic 
mechanisms can be evaluated by having them competing 
which each other Bifurcations appear when these two 
autocatalytic processes are placed in a nonequilibrium 
open-flow system, both being fed by the same incoming 
compound and with cross-inhibition between them: 



A; 



■*A 



A^B2 
Bi+B2^(P) 



(incoming flux) (31) 

(Direct AC) (32) 

(Autoinduced AC) (33) 

(cross inhibition) (34) 

(outgoing flux) (35) 

(outgoing flux) (36) 



In the case of total symmetry between and Bj, with 
the same direct autocatalystic mechanism, this system 
would correspond to the classical Frank model for the 
emergence of homochirality.^^ Because of the system sym- 
metry, the same probability to end up with either Bj or Bj 
is observed. 

The kinetic equivalence between template autocatalysis 
and autoinductive autocatalysis can be shown by mak- 
ing these two mechanisms to compete, replacing eq. ( [32) l 
and eq. ([33]) by the corresponding mechanism. Kinetic 
parameters have first been normalized so that each re- 
action leads on their own to the same kinetic behavior 
(sigmoidal evolution, half-reaction at 10^ s), and then 
multiplied by respectively a and P parameters in order 
to tune the respective velocity of each mechanism. The 
result is actually symmetrical between the two processes 
and only the fastest product is maintained in the system: 
Bj when a> j3, and Bj when a< [3 (see fig.[5](a)). As a 
consequence, while mechanistically different, these two 
autocatalysis are shown to be dynamically equivalent. 

This selectivity is independent of the relative stability 
of Bj and Bj, but is only possible for kinetics that are well 
adapted to the global influx of matter For slow kinetics. 
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5 Conclusion 




(a) Sharp bifurcation depending on the relative values of a and 
P for moderate reactivities. 




(b) Different zones of behaviors: majority of A for a, /3 <C 1^ ma- 
jority ofBi for a> P, majority of B2for a < P, and coexistence 
ofBj and 82 for a,p » 1. 

Figure 5 - Competition between template and autoinductive 
autocatalysis, generating respectively and B2 compounds 
from the same A compound. Incoming flux of A, and outgoing 
fluxes of and B2, 10"= M.s'^ = I, K^^ = K^^ = 
100. Direct autocatalysis: T^c = 10"^ ■ a, T^c = lO"** ■ a. 
Autoinduction, according to fig. ^d): Fj = p, r2 = = 



system and of self maintenance of the whole network by 
crossing energetical fluxes.E^E^l Confusion among these 
different phenomena can be pinpointed in the literaturejSl 
when the failure of autoinductive sets to be maintained 
do not originate from a difference of behavior between 
autocatalytic and autoinductive mechanisms, but from a 
defect in the closure of the system (e.g. induced by the 
leakage of some components). 



5 Conclusion 

Important distinctions need to be made between mech- 
anistic and dynamic aspects of autocatalysis. One single 
mechanism can produce different dynamics, while iden- 
tical dynamics can originate from different mechanisms. 
Thus, a pragmatic definition of autocatalysis have to be 
based on a kinetic signature, in order to classify the sys- 
tems according their observable behavior, rather than on 
a mechanistic signature, that would instead classify the 
systems according to the origin of their behavior. All the 
different autocatalytic processes described in this work 
are able to generate autocatalytic kinetics. They can con- 
stitute a pathway towards the onset of "self-sustaining au- 
tocatalytic sets", as chemical attractor in non-equilibrium 
networks. However, the problem of the evolvability of 
such systems must be kept in mind."*^ If a system evolves 
towards a stable attractor, no evolution turns out to be 
possible. There is the necessity of "open-ended" evolu- 
tion"*"* i.e. the possibility for a dynamic set to not only 
maintain itself (i.e. as a strict autocatalytic system) but 
also to act as a "general autocatalytic set", redounding 
upon the concept originally introduced by Muller® for the 
autocatalytic power linked to mutability of genes. For 
example, insights can be gained by a deeper and renewed 
study of the evolution of prions as a simple mechanism of 
mutable autocatalytic systems. 



100 • p, Kr 
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there is a flush of the system, and neither nor can 
be maintained. For fast kinetics, the system is close to 
equilibrium, the compounds Bj and Bj being both present 
in proportion to their respective stability (see fig.jsj^b)). 
Such result is well known for open flow Frank systems.'SSj 

4.3 From Autocatalytic Processes 
towards Autocatalytic Sets 

These competitive systems are able to dynamically 
maintain a set of components, to the detriment of others. 
These autocatalytic networks must however not be con- 
fused with autocatalytic sets. This latter notion is rather 
popular in the artificial life literature, but relies much 
more on the cooperation between autocatalytic mecha- 
nisms than on the competition that has just been detailed 
here. This implies a notion of material closure of the 



The kinetic behavior of three different mechanisms for 
autocatalytic transformations have been studied in details. 
The methodology consists in establishing the different 
chemical fluxes of the network. The relationship between 
these fluxes can be simplified by assuming the QSSA for 
relevant compounds. The purpose is then to establish the 
expression of the transformation flux (p as a function of 
the concentration of the reactants and the products. 

6.1 Indirect autocatalysis 

The four fluxes of fig.|3][c) can be written as: 

if, = r,(,V^Vo - Vc) (37) 

^2 = TaCVc - VsVe) (38) 

V3 = TsCVh - Vb) (39) 

^4 = ^4^Vb - Vd) (40) 
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The QSSA for D comes down to if if 

r,v,,Vo-r,Vc = r^Vs-r^Vo 
[r,v^ + r^)Vo = r^VB + r,Vc 



riV^ + r4 



(41) 
(42) 

(43) 



Replacing Vj, by eq. ( 43 1 in eq. ( 37]) gives 



rir4 

riV^ + r4 ^ ''riV^ + r4 



rir4 



The QSSA for E comes down to '4>2 — f3' 

r2Vc-r2VBVE = r^VE-r^Vs 



(44) 

(45) 
(46) 

(47) 
(48) 



Replacing Vg by eq. ( [48] ) in eq. ([38j) by Eq. gives: 



9z = T2\Vc- Vi 



2\ 



(49) 

(50) 
(51) 



At last, the QSSA for C comes down to (/Jj ^ 1/72 = 



Combining eq. (461 and Eq, eq. (51 1 gives: 



Vr 



with r' 



and r' 



r' + r' 
^ 1 ^ ^ 2 



rir4 



(52) 
(53) 
(54) 



Replacing by eq. ( 52 1 in eq. ([46]) gives: 



^ 1 ^ ^ 2 



^ 1 ^ 2 



(55) 

(56) 
(57) 



Replacing r[ and by their expression given in eq. ( 53 1 
and eq. ([54]) then gives: 



^ ■ 



VaVb-V^ 



± + ± + + Ys. 

r, r, r4 r. 



(58) 



6.2 Autoinductive autocatalysis 

The three fluxes of fig.]3][d) are: 

^1 = r,(,v^VE - Vc) 

'f2=r2(.Vc-VBVE) 
f3 = ^3(.VbVe, - Vc) 



(59) 
(60) 
(61) 



The QSSA for C comes down to (fi + i^^ — ¥'2j ^^d the 
QSSA for E comes down to (f-^ :^ if 2- This implies that 
1/73 <i if I, so that with ip3 = £ and, ipi = if, we obtain: 



(p2 = ^ + £^^ 



(62) 



In that context, eq. (61 ) gives: 



Vc = VbVe. - 



(63) 



Combining eq. ( 59 1, eq. ( 60 1 in eq. ( 62 1 then gives: 



r2Vc-r2VsVE = r^v^VE-r,Vc + e (64) 

(Fi + r2)Vc - e 



(65) 



Replacing by its value given in eq. (63 1 leads to: 



(ri + r2)VH,VB 



(66) 
(67) 



The flux of destruction of A can be computed by replacing 



in eq. ( 59 1 by eq. ( 67 ) (computing the flux of formation 



of B from eq. ( 60 ) would of course give the same result) : 



^ , (£^+£2)}^ , 
</' = rj V^— — — — — VbVe' I (68) 



(69) 
(70) 



riV^ + TjVB 



rir2% 



VbVa-V 



riVA + r2VB 
2 



riV^ + r2VB 

The law of conservation of £ compounds leads to: 

rir2vO(VBV^-Vg2) 



{T^V^ + T2VB\l + rcVc + rEVE) 



(71) 



with rc = Kc/Kg, and = K^/Kp . Assuming that E* is 
the much more stable than C and E, rc and <. 1, so 
that we finally obtairj^ 



riVA + r2VB 



WA-yn) 



(72) 



^Without the hypothesis of a large stabiUty of E* , not neglecting the 
and terms eventually leads to add Vg terms to the denominator, 
which will tend to destroy the autocatalytic effect. 
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6.3 Iwamura's model 

The five fluxes of fig.|3][e) are: 

'■P2 = ^2(yPAVx - VpAx) 
^3 = ^siYpAX ~ Vp\Ax) 

'■P4 = r4(Vp|^ - VpV^) 

^5 = ^si^P\AX^A - ^PA^AX) 



(73) 
(74) 
(75) 
(76) 
(77) 



The QSSA for P leads to Lp^ = ip^; for PA, it leads to 
^2 = ^1 + V'S; for PAX it leads to = ¥'2; for P\AX, it 
leads to ip^ = p^+p^. The fluxes can thus be decomposed 
into two elementary fluxes: 



(/3i = e 


(78) 




(79) 


if^ = ip + e 


(80) 




(81) 


^5 = ^ 


(82) 



if is the flux of the catalytic reaction, and e the flux of 
the non-catalytic reaction, so that e <. 1^. This would 
typically be characterized by -C Fg. 
1^2 = leads to: 



P\AX ■ 



^2 ^2 
T; ^PAX ~ TT^PA^X 
i 23 i 3 



(83) 



WithF23 = F2F3/(F2+F3) 

— ^5 leads to: 



P\AX 



^PAX + 



F, + F,V, 



VpaV, 



AX 



(84) 



5 "A 



Combining eq. (831 and eq. (841, eliminating Vp|^ leads 
to: 



PAX 



^fVAx + Vx + Y^V^V^ 



"r23^PA 



Combining eq. (85 1 with eq. (83 1 leads to: 

^5Vax+^23Vx, 



P\AX 



r5VA+F23 



-V, 



PA 



= leads to: 



(85) 



(86) 



(87) 



Combining eq. (871 and eq. (861 leads to: 

^ ^ F4(F5 + F23 Vx) + TiCFsV^ + F23) 
' (riV^ + F4V^)(F5V^ + F23) 

The flux of production of AX can be computed from 



PA 



(88) 



eq. ([74]), eq. ([75]) or eq. ([77]), which leads to: 

VaVx - Vax 



f ■ 



^5^23^PA 



(89) 



Combining eq. (891 and eq. (881 leads to: 



rs(^v^ + v^)+r23(^ + ^ 



(90) 



This can be simplified in an ideal case, assuming that the 
compound P is the mores stable compound among P, PA, 
PAX and P|AX, so that Vp Vp, and that the reactivities 
are so that F^ -C (F4,F5) <C F23 (i.e. assuming that 
reaction 1 is very slow, and that reactions 2 and 3 are very 
fast), which leads to: 



VaVax 



^AX 
Vy 



(91) 



Acknowledgment: This work was done within the scope 
of the European program COST "System Chemistry" CM0703. 
We additionally thank R. Pascal for useful discussions. 



References 

[1 

[3: 

[4 

[5: 
[6: 



[7 

[8: 
[9: 
[10: 

[11 

[12: 
[13: 

[14 

[is: 



23 



Ostwald, W Uber autokatalyse. Ber. Verh. Kgl. Sachs. Ges. 
Wiss. Leipzig, Math.- Phys. Classe 42, 189-191 (1890). 

Laidler, K. J. The development of theories of catalysis. 
Arch. Hist. Exact Sci. 35, 345-374 (1986). 

Ostwald, W. Lehrbuch der Allgemeinen Chemie, 263 (En- 
gelmann, Leipsic, 1902), 2nd edn. 

Ostwald, W. Outlines of general chemistry (trad. Taylor, 
W.W.), chap. XI.l, 301 (Macmillan and co., 1912). 

Lotka, A. J. Contribution to the theory of periodic reac- 
tions. J. Phys. Chem. 14, 271-274 (1910). 

Robertson, T. Further remarks on the normal rate of 
growth of an individual, and its biochemical significance. 
Devel. Genes Evol. 26, 108-118 (1908). 

Lotka, A. J. Analytical note on certain rhythmic relations 
in organic systems. Proc. Natl. Acad. Sci. USA 6, 410 
(1920). 

MuUer, H. J. Variation due to change in the individual 
gene. Amer. Naturalist 56, 32-50 (1922). 

Lotka, A. J. Elements of physical biology (Williams & 
Wilkins company, 1925). 

Witzemann, E. J. Mutation and adaptation as component 
parts of a universal principle: IL the autocatalysis curve. 
Amer Naturalist 67, 264-275 (1933). 

Malcai, O., Biham, O., Richmond, R & Solomon, S. Theo- 
retical analysis and simulations of the generalized Lotka- 
Volterra model. Phys. Rev. E 66, 031102 (2002). 

Frank, F. C. Spontaneous asymmetric synthesis. Biochem. 
Biophys. Acta 11, 459-463 (1953). 

Plasson, R., Kondepudi, D. K., Bersini, H., Commeyras, A. 
& Asakura, K. Emergence of homochirality in far-from- 
equilibrium systems: mechanisms and role in prebiotic 
chemistry. Chirality 19, 589-600 (2007). 

Kondepudi, D. K., Kaufman, R. J. & Singh, N. Chiral sym- 
metry breaking in sodium chlorate crystallization. Science 
250, 975-976 (1990). 

Soai, K., Shibata, T, Morioka, H. & Choji, K. Asymmetric 
autocatalysis and amplification of enantiomeric excess of 
a chiral molecule. Nature 378, 767-768 (1995). 



9 



References 



[16] Plasson, R. Comment on "re-examination of reversibil- 
ity in reaction models for the spontaneous emergence of 
homochirality". J. Phys. Chem. B 112, 9550-9552 (2008). 

[17] Blackmond, D. An examination of the role of autocatalytic 
cycles in the chemistry of proposed primordial reactions. 
Angew. Chem. 48, 386-390 (2009). 

[18] Peusner, L., Mikulecky, D. C, Bunow, B. & Caplan, S. R. A 

network thermodynamic approach to hill and king-altman 
reaction-diffusion kinetics. J. Chem. Phys. 83, 5559-5566 
(1985). 

[19] Mikulecky, D. C. Network thermodjmamics and complex- 
ity: a transition to relational systems theory. Comp. Chem. 
25, 369-391 (2001). 

[20] Plasson, R. & Bersini, H. Energetic and entropic analy- 
sis of mirror symmetry breaking processes in a recycled 
microreversible chemical system. J. Phys. Chem. B 113, 
3477-3490 (2009). 

[21] von Kiedrowski, G. Minimal replicator theory I: Parabolic 
versus exponential growth. Bioorg. Chem. Front. 3, 115- 
146 (1993). 

[22] Wills, R R., Kauffman, S. A., Stadler, B. M. R. & Stadler, B F. 
Selection djTiamics in autocatalytic systems: Templates 
replicating through binary ligation. Bull. Math. Biol. 60, 
1073 - 1098 (1998). 

[23] Mehra, A. et al. Circadian rhythmicity by autocatalysis. 
PLoS Comput. Biol. 2, e96 (2006). 

[24] Wang, Z.-X. & Wu, J.-W. Autophosphorylation kinetics of 
protein kinases. Biochem. J. 368, 947-952 (2002). 

[25] Iwamura, H. et al. Probing the active catalyst in product- 
accelerated proline-mediated reactions. J. Am. Chem. Soc. 
126, 16312-16313 (2004). 

[26] Ashkenazi, M. & Othmer, H. G. Spatial patterns in coupled 
biochemical oscillators. J. Math. Biol. 5, 305-350 (1977). 

[27] Nielsen, K., Sorensen, R G. & Hynne, F. Chaos in glycolysis. 
J. Theor. Biol. 186, 303 - 306 (1997). 

[28] Eigen, M. & Schuster, R The hypercycle, a principle of 

natural self-organization. Naturwiss. 64, 541-565 (1977). 

[29] Flores, R. et al. Viroids: the minimal non-coding RNAs 
with autonomous replication. FEBS Lett. 567, 42 - 48 
(2004). 

[30] Szathmary, E. Simple growth laws and selection conse- 
quences. Trends Ecol. Evol. 6, 366-370 (1991). 



[31] Nowak, M. A. Evolutionary Dynamics: exploring the equa- 
tions of life. (Harvard Univ. Press, 2006) . 

[32] von Kiedrowski, G. A self-replicating hexadeox3mucleo- 
tide. Angew. Chem. 25, 932-935 (1986). 

[33] Szathmary, E. & Gladkih, I. Sub-exponential growth and 
coexistence of non-enz5miatically replicating templates. J. 
Theor Biol. 138, 55 - 58 (1989). 

[34] Lifson, S. & Lifson, H. A model of prebiotic replication: 
Survival of the fittest versus extinction of the unfittest. J. 
Theor Biol. 199, 425 - 433 (1999). 

[35] Scheuring, I. 8;; Szathmary, E. Survival of replicators 
with parabolic growth tendency and exponential decay. J. 
Theor Biol. 212, 99 - 105 (2001). 

[36] Wagner, N. 8;; Ashkenasy, G. SjTnmetry and order in 
systems chemistry J. Chem. Phys. 130, 164907 (2009). 

[37] Gridnev, I., Serafimov, J., Quiney, H. 8;; Brown, J. Reflec- 
tions on spontaneous asymmetric sjmthesis by amplifying 
autocatalysis. Org. Biomol. Chem. 1, 3811-3819 (2003). 

[38] Schiaffino, L. 8;; Ercolani, G. Unraveling the mech- 
anism of the soai asymmetric autocatalytic reaction by 
first-principles calculations: Induction and amplification 
of chirality by self-assembly of hexamolecular complexes. 
Angew. Chem. 120, 6938-6941 (2008). 

[39] Cruz, J. M., Parmananda, R & Buhse, T. Noise-induced 
enantioselection in chiral autocatalysis. J. Phys. Chem. A 
112, 1673-1676 (2008). 

[40] Kauffman, S. A. Autocatalytic sets of proteins. J. Theor 
Biol. 119, 1-24 (1986). 

[41] Hordijk, M. & Steel, M. Detecting autocatalytic, self- 
sustaining sets in chemical reaction systems. J. Theor BioL 
227, 451-461 (2004). 

[42] Benko, G. et al A topological approach to chemical orga- 
nizations. Artificial Life 15, 71-88 (2009). 

[43] Vasas, V, Szathmary, E. & Santos, M. Lack of evolvabil- 
ity in self-sustaining autocatalytic networks constraints 
metabolism-first scenarios for the origin of life. Proc. NatL 
Acad. Sci. USA 107, 1470-1475 (2010). 

[44] Ruiz-Mirazo, K. Enabling conditions for 'open-ended 
evolution'. Biol. Philos. 23, 67 (2007). 

[45] Li, J., Browning, S., Mahal, S. R, Oelschlegel, A. M. & 
Weissmann, C. Darwinian evolution of prions in cell 
culture. Science 327, 869-872 (2010). 



10 



